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Abstract 

We present four continuations of the critical nonlinear Schrodinger equation (NLS) beyond the sin- 
gularity: 1) a sub-threshold power continuation, 2) a shrinking-hole continuation for ring-type so- 
lutions, 3) a vanishing nonlinear-damping continuation, and 4) a complex Ginzburg-Landau (CGL) 
continuation. Using asymptotic analysis, we explicitly calculate the limiting solutions beyond the 
singularity. These calculations show that for generic initial data that leads to a loglog collapse, 
the sub-threshold power limit is a Bourgain-Wang solution, both before and after the singularity, 
and the vanishing nonlinear-damping and CGL limits are a loglog solution before the singularity, 
and have an infinite- velocity expanding core after the singularity. Our results suggest that all NLS 
continuations share the universal feature that after the singularity time Tc, the phase of the singular 
core is only determined up to multiplication by e'^. As a result, interactions between post-collapse 
beams (filaments) become chaotic. We also show that when the continuation model leads to a point 
singularity and preserves the NLS invariance under the transformation t — )■ — t and ip ^ ip* , the 
singular core of the weak solution is symmetric with respect to Tc- Therefore, the sub-threshold 
power and the shrinking-hole continuations are symmetric with respect to Tc, but continuations 
which are based on perturbations of the NLS equation are generically asymmetric. 



1. Introduction 

The focusing nonlinear Schrodinger equation (NLS) 

#t(t,x) + AV' + |V|'"V' = 0, Vo(0,x) =^o(x) Gi?\ (1) 

where x = (xi, x^) G M"^ and A = dx^xi + • • 'dx^xa^ is one of the canonical nonlinear equations in 
physics, arising in various fields such as nonlinear optics, plasma physics, Bose-Einstein condensates 
(BEG), and surface waves. In the two-dimensional cubic case, this equation models the propagation 
of intense laser beams in a bulk Kerr medium. In that case, is the electric field envelope, t is 
the direction of propagation, d = 2, xi and X2 are the transverse coordinates, and a = 1 (cubic 
nonlinearity) . 

In 1965, Kelley showed that the two-dimensional cubic NLS admits solutions that collapse 
(become singular) at a finite time (distance) Tc [1]. Since physical quantities do not become 
singular, this implies that some of the terms that were neglected in the derivation of the NLS, 
become important near the singularity. Therefore, the standard approach for continuing the solution 
beyond the singularity has been to consider a more comprehensive model, in which the collapse is 
arrested. 
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In this study, we adopt a different approach, and ask whether singular NLS sohitions can be 
continued beyond the singularity, within the NLS model. By this we mean that the solution satisfies 
the NLS both before and after the singularity, and a matching ("jump") condition at the singularity. 
The motivation for this approach comes from hyperbolic conservation laws, where in the absence 
of viscosity, the solution can become singular (develop shocks). In that case, there is a huge 
body of literature on how to continue the inviscid solution beyond the singularity, which consists 
of Riemann problems, vanishing-viscosity solutions, entropy conditions, Rankine-Hogoniot jump 
conditions, specialized numerical methods, etc. In contrast, two studies from 1992 by Merle [2, 3], 
and a recent study by Merle, Raphael and Szeftel [4], addressed this question in the NLS. Tao [5] 
proved the global existence and uniqueness in the semi Strichartz class for solutions of the critical 
NLS. Intuitively, these solutions are formed by solving the equation in the Strichartz class whenever 
possible, and deleting any power that escapes to spatial or frequency infinity when the solution 
leaves the Strichartz class. These solutions, however, do not depend continuously on the initial 
conditions, and are thus not a well-posed class of solutions. Stinis [6] studied numerically the 
continuation of singular NLS solutions using the t-model approach. 

In [2], Merle presented an explicit continuation of a singular NLS solution beyond the singularity, 
which is based on reducing the power (L^ norm) of the initial condition of the explicit blowup 
solution V'cxpiicit.a) scc (8), below the critical power for collapse Per- This continuation has two key 
properties: 

1. Property 1: The solution is symmetric with respect to the singularity time Tc- 

2. Property 2: After the singularity, the solution can only be determined up to multiplication 
by a constant phase term. 

Merle's breakthrough result, however, applies only to the critical NLS {ad = 2), and only to 
the explicit blowup solutions V'cxpUcit.a- Recently, Merle, Raphael and Szeftel [4] generalized this 
result to Bourgain-Wang singular solutions [7], i.e., solutions that have a singular component that 
collapses as '(/'explicit, a i and a non-zero regular component that vanishes at the singularity point and 
does not participate in the collapse process. 

In [3], Merle presented a different continuation, which is based on the addition of nonlinear 
saturation. This study showed that, generically, as the nonlinear saturation parameter goes to zero, 
the limiting solution can be decomposed beyond Tc into two components: A J-function singular 
component with power m{t) > Pa-, and a regular component elsewhere. Similar results follow from 
the asymptotic analysis of Malkin [8], which suggests that m(t) = Per for Tc < t < oo. 

It is thus useful to distinguish between two types of continuations: 

1. A point singularity, in which the weak solution is singular at Tc, but regular (i.e., in H^) 
for t > Tc. 

2. A filament singularity, in which the weak solution has a 5-function singularity for Tc < t < 
Tq, where Tc < Tq < oo. The motivation for this terminology comes from nonlinear optics, 
where collapsing beams can form long and narrow filaments (which, to leading order, can be 
viewed as an "extended" (5- function) . 

In this work we propose four novel continuations that lead to a point singularity, and obtain 
explicit formulae for the solution beyond the singularity. Our main findings are: 

1. The non-uniqueness of the phase beyond the singularity (Property 2) is a universal feature of 
NLS continuations. 

2. The symmetry with respect to the singularity time (Property 1) holds only if the continuation 
is time reversible and leads to a point singularity. Therefore, it is non-generic. 
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The paper is organized as follows. In section 2 we provide a short review of NLS theory. In 
section 3 we present Merle's continuation of V'cxpUcit.cn and illustrate it numerically. In section 4 we 
generalize this approach and present a sub-threshold power continuation, which can be applied to 
generic initial condition of the critical NLS. We compute asymptotically the limiting solution, and 
show that it is a Bourgain-Wang solution, both before and after the singularity. In particular, this 
continuation preserves the two key properties of Merle's continuation. In section 5 we show that 
Property 1 (symmetry with respect to Tc) holds for time-reversible continuations that lead to a 
point singularity. In Section 6 we discuss the nonlinear-saturation continuation, which leads to a fil- 
ament singularity. In section 7 we show that because of the phase non-uniqueness (Property 2), the 
interaction between post-collapse beams is chaotic. In section 8 we present a vanishing-hole contin- 
uation, which is suitable for ring-type singular solutions. This continuation is time-reversible, and 
it satisfies Properties 1 and 2. In section 9 we present a vanishing nonlinear-damping continuation, 
and compute the continuation asymptotically in two cases: 

1. The nonlinear-damping continuation of the explicit solution V'expiicit.a is, up to an undeter- 
mined constant phase, given by V'expiicit.Ka with k 1.614. 

2. The nonlinear-damping continuation of solutions that undergo a loglog collapse has an infinite- 
velocity expanding core, with an undetermined constant phase. 

Therefore, the phase becomes non-unique beyond the singularity (Property 2). In contrast with 
previous continuations, however, the solution is asymmetric with respect to the singularity time 
(i.e., the solution does not satisfy Property 1). This is to be expected, as the nonlinear-damping 
continuation is not time-reversible. In Section 10 we show that the continuation which is based on 
the complex Ginzburg-Landau (CGL) limit of the NLS, is equivalent to the vanishing nonlinear- 
damping continuation. In section 11 we show a continuation of singular solutions of the linear 
Schrodinger equation. In this case, the limiting phase beyond the singularity is unique. This 
shows that the post-collapse non-uniqueness of the phase (Property 2) is a nonlinear phenomenon. 
Section 12 concludes with a discussion. 

1.1. Level of rigor 

The results which are derived in this manuscript are non-rigorous, and are based on asymptotic 
analysis, numerical simulations, and physical arguments. To emphasize this, we use the terminology 
Continuation Results, rather than Propositions or Theorems. 

2. Review of NLS theory 

We briefly review NLS theory, for more information see [9, 10, 11]. The NLS (1) has two 
important conservation laws: Power conservation^ 





and Hamiltonian conservation 






We call the norm the power, since in optics it corresponds to the beam's power. 
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The NLS admits the waveguide solutions ip = e^^R{r), where r = |x|, and R is the solution of 

R"{r) + ^^R' - R + R^''+^ = 0, R'{0) = 0, i?(oo) = 0. (3) 
r 

When d = 1, the solution of (3) is unique, and is given by 

R^ix) = (1 + a)^/'^''cosh~^^''{ax). (4) 

When d > 2, equation (3) admits an infinite number of solutions. The solution with the minimal 
power, which we denote by is unique, and is called the ground state. 

When ad < 2, the NLS is called subcritical. In that case, all solutions exist globally. 
In contrast, both the critical NLS {ad = 2) and the supercritical NLS {ad > 2) admit singular 
solutions. 

Let ip{t,x.) be a solution of the NLS (1). Then, remains a solution of the NLS (1) under the 
following transformations: 

1. Spatial translations: "0(^5 — ?• ip{t,:K + xq), where xq G M*^. 

2. Temporal translations: ^{t,x) — t- ■ip{t + to,x.), where to G 

3. Phase change: il^{t,:x.) — )• e'^^ip{t,:x.), where € M. 

4. Dilation: ^A(t,x) A^/'^V(A^t, Ax), where A G IR+. 

5. Galilean transformation: 0(t,x) — ip{t,:K — ct)e*'^'''/^~*l'^l^*/^, where c G M.'^. 

Therefore, multiplying the initial condition by a constant phase e*^ does not affect the solution. 
In addition, by the Galilean transformation, multiplying the initial condition by a linear phase 
term 'ipo{x) — ?• V'o(x)e*^'^/^ does not affect the dynamics, but rather causes the solution to be tilted 
in the direction of n = (1, c) G M x M''. 

2.1. Critical NLS 

In the critical case ad = 2, equation (1) can be rewritten as 

iVt(t,x) + A0 + |0|4/'^0 = O, Vo(0,x) =Vo(x) G/?\ (5) 

and equation (3) can be rewritten as 

R"{r) + ^^R' - R + R^/'^+^ = 0, R'{0) = 0, i?(oo) = 0. (6) 
r 

Theorem 1 (Weinstein [12]). A sufficient condition for global existence in the critical NLS (5) is 
that 11^0 111 < -fcr; where Per = \\R^^M\2! ^^^^ -^^^^ the ground state of equation (6). 

The critical NLS (5) admits the explicit solution 

r) = ^fi<»'(^) (7a) 

where 



L{t) =T,-t, T{t) = I ds = (7b) 
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More generally, applying the dilation transformation with \ = a and the temporal translation 
Tc — 7- a^Tc shows that the critical NLS (5) admits the explicit solutions 



V'explicit,« (^5 ^) — ^/r, ^ 



1 „^n^ f r \ I L. 
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where 

/■* 1 11 

= a{Tc - t), Ta{t) = / —y-- ds = -, a > 0. (8b) 



Even more generally, by the Galilean transformation, 

V'eipLt.Ji, X) = Vexplicit,.(t, X - C • t)e--/2-l-l^*/4^ (9) 

are also explicit solutions of the critical NLS (5). 

The explicit solutions (7)-(9) become singular at t = Tc. These solutions are unstable, however, 
as they have exactly the critical power for collapse. Therefore, any infinitesimal perturbation which 
decreases their power, will arrest the collapse. 

The only minimal-power blowup solutions (i.e., singular solutions whose power is exactly Per) 
are given by ^^'^pl^j^.a^ 

Theorem 2 (Merle [2, 3]). Let if) he a solution of the critical NLS (5) which blows up at a finite 
time Tc > 0, such that HV'olli — -^cr- Then, there exist a G M"*", G M, and xo,c G M.'^, such that 
forO<t< Tc, 

V'(t,x)=<'*',^,,„(t,x-xo)e^^ (10) 

When an NLS solution whose power is slightly above Per undergoes a stable collapse, it splits 
into two components: A collapsing core that approaches the universal V'ijCo) profile and blows up at 
the loglog law rate, and a non-collapsing tail {(f)) that does not participate in the collapse process: 

Theorem 3 (Merle and Raphael [13], [14], [15], [16], [17], [18], [19]). Let d = 1,2,3,4,5, and let ^ 
he a solution of the critical NLS (5) that hecomes singular at Tc- Then, there exists a universal 
constant m* > 0, which depends only on the dimension, such that for any ipo G LL^ such that 



Im J ip^ViJo\ 
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^'cr< llV'olli <^'cr + m*, Hci^Po) := g(V'o) " ( "'^i.:,,' ) < 0, 

the following hold: 



1. There exist parameters (T(t), xo(i), L(t)) G M x M*^ x M"*", and a function 7^ </> G L"^, such 
that 



where 



1^2 

ipit, x) - V'k(o) {t, X - xo(t)) — > 0(x), t — > Tc, 



and R^^^ is the ground state of equation (6). 
2. Ast — > Tc, 



^'"~^( iog|iS(r!-t)| )"' CosX'^i-)- (11) 
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NLS solutions whose power is slightly above Per can also undergo a different type of collapse, 
in which the collapsing core approaches the explicit blowup solution V'cxpUcit.a and blows up at a 
linear rate, but the solution also has a nontrivial tail {(j)) that does not participate in the collapse 
process: 

Theorem 4 (Bourgain and Wang [7]). Let d = 2, let Aq he a given integer, and let A > Aq he a 
large enough integer. Let 4> £ Xa = {f £ with (1 + \x\^)f £ L'^}, and let z £ C{{T* ,Tc\, X^) 
he the solution to the critical NLS (5), subject to z{t = Tc) = (j), where T* < Tc is the maximal 
time of existence of z. Assume that (p vanishes to high order at the origin, i.e., D'^(j){0) = for 
\a\ < A — 1. Then, there exists T* <to < Tc and a unique solution tpBW G C{[to,Tc), Xaq) to (5), 
such that 

IIV'Bw(t) - V'explicit,a(t) " ^(t)lk^, < \Tc " (12) 

The Bourgain- Wang solutions ipEwii, x) are unstable [4], since their singular core is the unstable 

blowup solution V'cxplicit,a- 



3. Merle's first continuation 

In [2], Merle presented a continuation of the explicit blowup solution V'expUdt.a beyond the 
singularity. To do that, he considered the solution x/j^^^t, r) of the critical NLS (5) with the initial 
condition 

4'^(r) = (1 - e)V'cxpiicit,a(i = 0,r), < e < 1. (13) 

Since the power of ■0^^-' is below Per, it exists globally. Therefore, it is possible to continue the 
singular solution 0expiicit,o beyond the singularity, by taking the limit of 0^^-* as e — )■ 0+. The 
limiting solution is given in the following Theorem: 

Theorem 5. [2] Let ip^^\t,r) he the solution of the critical NLS (5) with the initial condition (13). 
Then, for any G R, there exists a sequence e„ — )■ 0+ ( depending on 9 ), such that 

where 'i/'cxpiicit(i, t) is given hy (8). 

Theorem 5 shows that after the singularity, the limiting solution is completely determined, up 
to multiplication by e*^, and is symmetric with respect to Tc, i.e., 

lim '(/'(^")(rc + t,r) = e*^ lim 0*(^")(rc - 1, r), t > 0. (15) 

£n^0+ e„-s>0+ 

3.1. Simulations 

In order to provide a numerical illustration of Theorem 5, let il)^^\t,x) be the solution of the 
one-dimensional critical NLS 

i0f(t,x) + 0^.^. + |0|V = O, (16a) 

with the initial condition 

4'^(x) = (l-e)0expiicit(O,x) = (l-e)^p(^^^ e'^-^^, Te = 0.25. (16b) 
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Let us define the width of the solution ^l)^^\t,x) as 



Le{t) :-- 



RiO) 



V'(<^)(t,o) 



2/d 



(17) 



see (7a) and (25). Figure 1(a) shows that lim£_j>o+ Li;{t) = \Tc — t\, both for t < Tc and for Tc < t, 
in agreement with Theorem 5. 




arg v<^'{t,0). 




(b) 



0.5 



1 1.5 

t/T 



2.5 



Figure 1: Solution of (16). (a): L,{t) for e = 5 • lO"*, 2.5 • 10"* and lO"*. Solid line is L = \Tc - 1\. (b): Accumulated 
phase at a; = for e = 0.9 • 10"^, 10"^, and 1.1 ■ 10"^. Solid line is arg(V'expiicit(t, 0)). 



In order to observe the loss of the phase after the singularity, we compute the effect of small 
changes in the initial condition on the phase, by solving equation (16) with e = 0.9 • 10~^, 10~^ 
and 1.1 • 10~^. Figure 1(b) shows that these 0(10"^) changes in the initial condition lead to 0(1) 
changes in the phase for t > Tc, which is a manifestation of the post-collapse phase loss as e — )• 0+. 

4. Sub-threshold power continuation 

The main wcakiiGSS of ThGorcm 5 is that it only applies to the explicit solutions V^expiicitja* 

We 

now generalize Merle's continuation to generic initial profile F{x) G H^, as follows. Consider the 
solution tjj^t, x; K) of the critical NLS (5) with the initial condition 

'0o(x; K) = K ■ F(x), F(x) e H\ K > 0. (18) 

Let 

= in{{K I ip{t,x;K) collapses at some < Tc{K) < oo}. (19) 

Let us consider the case where the infimum is attained, i.e., when 'ip{t;x; K^fi) becomes singular at 
a finite time. By construction, the solution of the critical NLS (5) with the initial condition 

= (1 - e) • ^^^(x), ^^^(x) = Kt, . F(x), (20) 

exists globally for < e ^ 1, but collapses for — 1 ^ e < 0. Therefore, as in Theorem 5, we 
can define the continuation of ip{t,x; Kth) beyond the singularity, by considering the limit of x/j^^^ 
as e — 7- 0+. Using asymptotic analysis, in Section 4.2 we derive the following result, which is the 
non-rigorous asymptotic analog of Theorem 5: 
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Continuation Result 1. Let il)^''\t,r) he the solution of the critical NLS (5) with the initial 
condition (20), where F{r) is radial. Assume that il)^'^^^\t,r) becomes singular at < Tc < oo. 
Then, for any 9 G M., there exists a sequence £n — ^ 0-|- (depending on ) and a function (p L , 
such that 

limt^T,- [lim£^o+ V'^^H^i - V'cxplicit,a(i, r)e*''o] = 0(r) 

[lim,„^o+ V'^^") {t, r) - VcVicit, J2r, - t, r)e*^T ^^^^ 



lim 



where the above limits are in L^, ^explicit, o is given by (8), Tc, a G and 6q € M. Therefore, 
locally near the singularity, the limiting solution satisfies Properties 1 and 2. 
In particular, the limiting width of the collapsing core is given by 

:= lim L(t;e) ~ alt - TJ, t^Tc±. (22) 

e-s>0+ 

By Theorem 2, if -(/'(t; x; Kth) is a singular solution which is not given by (10), then | ['(/'(t; x; ETj/j)] I2 > 
Per- In that case, 

U\\l = hm ||V^(^)||2-P„ = ||^(t;x;KtO||2-P„>0. 
e— >-0+ 



CoroUary 1. // 



the limiting solution in Continuation Result 1 is a Bourgain-Wang solution ipBW: both before and 
after the singularity. 

4.1. Simulations 

In order to illustrate the results of Continuation Result 1, we solve numerically the one- 
dimensional critical NLS (16a) with the initial condition 

4\x) = {l-e)- 4'"\x), ^'^(x) = KtHe--\ (23) 

We first compute the value of Kf^- We solve the NLS (16a) with the initial condition tpQ = Ke~^ . 
Figure 2(a) shows that for K = 1.48140 the beam collapses, while for K = 1.48139 the collapse is 
arrested. Therefore, Kth = 1.481395 ± 5 • 10"^. 

In Figure 2(b) we plot the solution of (16a) with the initial condition (23), as a function of t. The 
solution amplitude increases up to t = T^^^ := argmax^ \\'ilj^'^\t,x = 0)||oo, and then decreases. As 
expected, both T^^x and the maximal amplitude \'ip^'^\T^^^,0)\ increase as e — )■ 0+. In Figure 2(c) 
we plot T^^^ as a function of e. Extrapolation of these values to e = shows that 

Tc:= lim r^,, p. 8.00. (24) 

Since Tc < 00, ilj{t,x; Kth) is a singular solution. 

In order to confirm that the profile of the collapsing core is a rescaled R^'^^ profile, see equa- 
tion (21), in Figure 2(d) we plot l^'*-'^-'! for £ = 2 • 10~^ at t = 10 ^ 1.35T^^^, i.e., after it collapse 
has been arrested, and observe that for < x/L{t) < 5, the rescaled profile is indistinguishable 
from i?^''), while for 6 < x/L(t) the two curves are different. This confirms that the inner core 
collapses with the ^/'expiicit,a profile, but the outer tail does not. 

In Figure 2(e) we plot the solution width -^(.(t) for various values of e. The extrapolation of 
these curves to e = is in good agreement with the predicted linear limit (22) for 0.25 < t/Tc < 1.2, 
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Ql . 1 1 I . 1 I . _ 

5 5 10 12 




Figure 2: Solution of tlie one-dimensional critical NLS (16a) with: (a) The initial condition (18) with K = 1.48139 
(solid) and K — 1.48140 (dashes), (b)-(g) the initial condition (23) with Kth ~ 1.481395 with various values 
of s. (b) Solution for e = 2 • 10"^4 ■ 10"^ and 6 • 10"^ (c) The maximum focusing time T^^^. (d) Rescaled 
solution for e — 2 ■ 10~^ at t = 10 « l.S5T^ax (solid), and the R{x) profile (dashes), on a semi-logarithmic scale, 
(e) L^{t) for e = 6 • 10~^ 4 • 10"^ 2 • 10"'^ and 1 ■ 10"^ (f) Extrapolation of {L^(t)} from (e) to e = (solid). Dashed 
line is L = a\Tc — t\ with a — 0.0777 and Tc — 7.96. (g) Accumulated phase at a; = 0. 



see Figure 2(f). Finally, in order to observe the loss of the phase after the singularity, we compute 
the effect of 0(10"^) changes in the initial condition on the phase, by solving the one-dimensional 
NLS (16a) with the initial condition (23) with e = 2-10^^, 4-10"^ and 6-10"^. Figure 2(g) shows that 
these 0(10"^) chan ges in the initial condition lead to 0(1) changes in the phase for 8 < t < 10. Since 
Tc w 8, see (24), these 0(1) changes occur after the singularity, in agreement with Continuation 
Result 1. 

4-2. Proof of Continuation Result 1 
4-2.1. Adiabatic collapse - review 

In order to compute asymptotically the limit of V''' as e — )• 0-|-, we recall that, in general, the 
collapse of radial solutions with power close to Per can be divided into two stages, see [9]: 

1. During the initial non-adiabatic self- focusing stage, the solution "splits" into a collapsing 
core ^/'core and a non-collapsing "tail", i.e., 

where L{t) is the collapsing core "width" and po = 0(1). 
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core approaches the self-similar profile 

-/-H... = S = rW + ^!^. r = ^'^&, (25b) 

where R^^^ is the ground state of (6). In contrast, the "tail" continues to propagate forward. 
In particular, limj^j'^ ^/^taii = fpi^) ^ -^^i see Theorem 3. 

Once the profile of ipcom is close enough to ipRW , the dynamics of the collapsing core becomes 
nearly adiabatic, and is governed, to leading order, by the reduced equations [20, 21, 22] 

m = Lu = (26) 



where 



and 



2A| 

M '■ 



1 foa 

= lim eV('^-i)/2^(0)(r), M=- r'^{R^^^fr'^-^dr. (27b) 



In addition, the parameter (3 is proportional to the excess power above Per of the collapsing 
core -^core, i-e., 

o ||V'core||2 ~ -^cr 

/? , (28) 

see [8, 9]. 

4- 2. 2. Asymptotic analysis 

For simplicity, we assume that the initial condition (20) is radial. By construction, for the initial 
condition (20), 

^core 1 1 2 ^ ^cr £ < , 
^core||2 ^ -^cr £ > 0. 

Therefore, by (28), 

' /3 > £ < 0, 



/3 < £ > 0. ^^^^ 

Hence, by (27a), v{l3) = for £ > 0. Therefore /3t{t;e) = 0. Thus, the self-focusing dynamics is 
governed by 

Lu{t;e) = -^, (30) 

where /?(£) is independent of t. 

Let tad(^) > denote the time at which ^("^^ "enters" the adiabatic stage, i.e., when ip^^^ ~ '4'r(o) , 
so that equations (26) hold. Therefore, for t > tadi^), the dynamics is given by 

Luit) = -^, L{tad) = Lad{e), Ltitad) = LUe)- (31) 
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Lemma 1. Let L{t;e) he the solution of (31). Denote 

tad:= lim tadie), a := lim Lt{tad{ey, e), T, := lim ^-die)L'^aie) ^ 



Then, 

Proof. Two integrations show that the solution of (31) is given by 



Hm L(t; e) = alt — TJ, tad^t<oo. (33) 

£-5>0+ 



V ci(e) J ci(e) 



Lldi^) ' 

By (29), Hme^o+/3(e) = 0. Therefore, Um£^o+ci(e) = a^. Hence, Hme^o+ ^^(*; e) = a^{t-Tcf. □ 



Corollary 2. Um£_s.o+ ''"(^ = Tc]£) = oo. 
Proof. By (25b) and (33), 

f^" dt f'^" dt 
hm r(t = Tc\e) = hm / —77- — - = / -77 = cxo. 

e^o+ ^ ' s^o+Jq L^{t-e) Jo a^{T,-tY 

□ 

In order to go back from L[t) to ip, let us note that: 
Lemma 2. Let ^ji(o) be given by (25b). If L{t) = a\Tc — t\, then 

^«-(*'^) = ie-^^,V,,,„(2T.-t,r), t > T,, ^3^) 

where il>cxphcit,a{t,^) is given by (8), and 9q,6i E M. 

4-2.3. Proof of Continuation Result 1 

In Lemma 1 we saw that the solution L{t) of the reduced system is given by (33) for tad < t < 00. 
Therefore, by Lemma 2, when tad < t < Tc, 'ipR(o){t,r) = ipcxplicit,a{t,r)e^^° , and when Tc < 
t, V'R(o) it, r) = V':.piicit,a(2rc - r) lim,„^o e*'^(^") . 

Since arg^(t, 0) ~ arg ■0^{o) (t, 0) = T(t), Corollary 2 shows that the limiting phase becomes 
infinite at Tc, hence also for t > Tc- Therefore, for a given t > Tc and 9 £ M., there exists 
a sequence e„ — )• 0+, such that lim^^_^Q^ aigip^^"-\t,0) = 6. Since as t — )• Tc, V'core — V'ijco) 
and 'i/'tail ~^ <^(x) G L'^, the Proposition follows. 

5. Time-reversible continuations 

The continuation in Continuation Result 1 preserves Properties 1 and 2 of Merle's first continu- 
ation. We now show that these two properties hold for continuations of the NLS that preserve the 
NLS invariance under the transformation, 

t^-t and i^^ip*, (35) 
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and also satisfy some additional conditions. 

Continuation Result 2. Let 'ilj{t,x.) be a solution of the NLS (1) that blows up at T^, and 
let ip^^\t,x) be a smooth continuation of %l}(t^:si), such that 

1. ^/;('') exists globally for < e ^ 1. 

2. lim£^o+ V'^'''' (i, x) = x) in L'^''+'^ for < t < Tc. 

3. ip^^^ is invariant under the transformation (35). 

4. linit^Tc argV'(t,0) = oo. 

5. LetT^^^ := argmax^ r)||2(T+2 denote the time at which the collapse of^^^^ is arrested. 
Then, for all x G R'^, 

argV'(^)(r^,,,x) = a(^), a' e R. (36) 
Then, for any 6 €zR, there exists CI scQVLCTicG Sn — ^ O-j- (dcpcTidifiQ Oft 6 J , such that 

lim ^(^")(Tc + t,x) = e^V*(rc-t,x), i > 0. 

£n-5>0+ 

Hence, the continuation satisfies Properties 1 and 2. 

Proof. Assume first that a^^-* = 0. Then, from the invariance of ^^^^ under (35) it follows that 

V^(^) (T^,, + t, x) = r^^) (T^^^ - t, x), t > 0. (37) 
If a(^) / 0, then (37) holds for e'^^^'V^^^ i-e., 

e--'^V(^^(r^a. + t,x) = e-'^V*(^)(r;nax - t,x), t>0. (38) 

Hence, 

V'^'^(r^ax +i,x) = e*^*^V*^'^(r^ax - t,x), t > 0, (39) 

where 6^^^ = 2a^^\ Since lim£_^o+ = li™e-!>o+ arg^'^'^H^max; 0) = c», there exists Sn — 0+ 
such that 6 = lim£^_^o+ 0^^"^\ In addition, from limt_5.Tc ||V'(*)I|2(t+2 = oo and lime_5>o ||V'^(0ll2o-+2 = 
||V'(*)l|2(T+2 for < t < Tc, it follows that lim^^Q^ T^^^ = Tc. The result follows by taking the 
limit of (39). □ 

Continuation Result 2 holds for both the critical and the supercritical NLS. The interpretation 
of the conditions of Continuation Result 2 is as follows. Conditions 1 and 2 say that ip^^^ is a 
continuation of ip. Condition 3 says that the continuation is time reversible. Condition 4 says that 
the phase of the singular solution becomes infinite at the singularity. This condition holds for all 
known singular solutions of the critical and supercritical NLS. Condition 5 says that at i = T^^^, 
the solution is collimated. Intuitively, this is because the solution is focusing for t < T^^^^, and 
defocusing for t > T^^^. 

We now confirm that the sub-threshold power continuation of Continuation Result 1 satisfies 
the conditions of Continuation Result 2. By (25b), 

1 / r\ ,v +iMtrl dr^ 1 



Since L^{t) attains its minimum at T^^^, then L'^{T^^^) = 0. Therefore, avg ip^^\T^^^,r) = a^'^\ 
where a^^^ = t^{T^^^). Furthermore, since Ll{t) ~ ~ ||y^;jp; and since in the critical 
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NLS \\V^P\\l>M{T,-t)-^ [23], 



lim r^„^. = lim / j2>M = oo. 

£^0+ e^0+ L| Jq Tc-t 

Therefore, lim£__i.o+ ck^^-' = oo. Hence, lim£_5.o+ 6''''^-* = oo. 

Prom all the conditions of Continuation Result 2, the only one whose validity is questionable is 
condition 5. It is reasonable to expect that this condition would hold for the collapsing core. There 
is no reason, however, why it should hold for the non-collapsing tail. 

An immediate consequence of Continuation Result 2 is that 

Corollary 3. Under the conditions of Continuation Result 2, the limiting solution lim£j^_>.o+ ^/^(''") 
is in H^for t > T^. Hence, the continuation leads to a point singularity, and not to a filament 
singularity. 

In Section 6 we will see that time-reversible continuations can also lead to a filament singularity. 
In that case, however, condition 5 does not hold. 



6. Vanishing nonlinear-saturation continuation 

6.1. Merle's second continuation 

In [3] , Merle presented a different continuation, which is based on arresting the collapse with an 
addition of nonlinear saturation. 

Theorem 6. [3] Let d > 2, and consider radial initial data ipoii") £ fll^V'o £ -^^}; such that 
the solution 'ip{t,x.) of the critical NLS (5) blows up in finite time T^. For e > and 1 +A/d < q < 
{d + 2)/{d — 2), let ip£{t, x) be the solution of the saturated critical NLS 

iMt, x) + A-^ + IV'I^/'^V' - elV'l^" V = 0, ^^(0, x) = Mr)- (40) 

If for Tq > Tc, there is a constant C > such that f |xp|'0£(To, x)pdx < C, then: 

1. There is a function tp{t, x) defined for t < Tq, such that for all rQ > 0, ip £ C([0, Tq), L^(|x| > 
ro)), and V'e(*,x) ■0(i,x) in C{[0,To), L'^{\x\ > ro)) as e ^ 0. 

2. For t < Tq, there is m{t) > such that |'!/'e(t,x)p — ^ m{t)6{x.) + |'i/5(t,x)p as e — ^ in the 
distribution sense. Furthermore, 

(a) If m{t) / 0, then ||'i/'£(i,x)||j|^i — +oo as e — and m{t) > Per. 

(b) Ifm{t) = 0, there is a constant c > such that for alls, ||V'£(t,x)||j|^i < c, and ipsit,x) 
'4>{t, x) in L^. 

3. For allt <To, m{t) + f \ip{t,:xi)\'^dx = /|V'o(x)pdx. 

Theorem 6 shows that the vanishing nonlinear saturation continuation can lead to a filament 
singularity. The condition J |xp|'0£(ro, x) dx < C is believed to hold generically. 

6.2. Malkins analysis 

In [8] , Malkin analyzed asymptotically the solutions of the saturated critical NLS (40) with d = 2 
and q = b. Malkin showed that initially, the solution follows the non-saturated NLS solution 
and self-focuses. Then, the collapse is arrested by the nonlinear saturation, leading to focusing- 
defocusing oscillations. During each oscillation, the collapsing core loses (radiates) some power. As 
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a result, the magnitude of the oscillations decreases, so that ultimately, the solution approaches a 
standing- wave solution of the saturated NLS. 

If we fix the initial condition i{jq and let e — > 0+, then as e decreases, the collapse is arrested 
at a later stage. For example, in Figure 3 we plot the solution of the saturated critical NLS (40) 
with d = 2 and q = 5 with the initial condition ipo{r) = 3.079e~^ , and observe that as e decreases, 
the collapse is arrested at a later stage, and the oscillations occurs at higher amplitudes. In addition, 
we observe that as t increases, the oscillations decreases. Hence, it is reasonable to assume that the 
amplitude of the limiting standing-wave increases as e decreases, and goes to infinity as e — > 0+. 
In addition, as e — > 0+, the power of the standing-wave of (40) approaches Per- Therefore, 
Malkin's analysis suggests that 

m{t) = Per, Tc < t < CO, 

i.e., after the singularity the limiting solution consists of a semi-infinite filament with power Per; 
and a regular part with power HV'olli ~ -fcr- 
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Ql . . , 1 

2 4 6 8 

t 



Figure 3: Solution of the saturated NLS (40) with d=2 and g = 5 for e = 0.5-10"^ (solid), and e = 2-10^^ (dashes). 
Here, tpoir) = 3.079e-'''. 



6.3. Importance of power radiation 

The above results of Merle and Malkin strongly suggest that the continuation of singular NLS 
solutions with a vanishing nonlinear-saturation generically leads to a filament singularity. Since the 
NLS with a nonlinear saturation is time reversible, these results seem to be in contradiction with 
Continuation Result 2, see Corollary 3. Note, however, that in Continuation Result 2 we assumed 
that the solution phase is constant at the time T^^^ where its collapse is arrested (Condition 5). If 
this condition were to hold for the solution of the saturated NLS, then by Continuation Result 2, 
the solution at t = 2T^^^ would be given by tpQ. As Figure 3 shows, this is not the case. 

The mechanism which enables the filamentation is the loss (radiation) of power from the collaps- 
ing core to the surrounding background. Indeed, in the absence of radiation, the collapsing core of 
the solution of saturated NLS undergoes periodic oscillations, rather than approaches a standing 
wave [8]. We stress that the constant phase condition does hold asymptotically for the collapsing 
core [8]. It does not, however, hold for the regular part of the solution (the "tail").^ 



^In the terminology of Theorem 3, it holds for '0^(0) , but not for 0. 
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7. Chaotic interactions 



In Sections 3-5 we saw that time-reversible continuations have the property whereby the phase 
becomes non-unique after the blowup time. A-priori, this phase loss should have no effect, since 
multiplying the NLS solution by e*^ does not affect the dynamics. Nevertheless, we now show that 
this phase loss can affect the interaction between two post-collapse beams (filaments). 

Consider first the initial condition 

V'o(x) = V'exgitlO' X T xo) = Vcxpiicit(0, X ^ xo)e±^^-(™)/2. (41) 
By (9), the solution of the critical NLS (5) with the initial condition (41) is given by 

Cp^cit(*> X T Xo) = ^explicit (t, X T Xo ± C • t)e±--(xTxo)/2-^|cPt/4_ 

Therefore, V'expiicit(^' explicit blowup solution (7), centered initially at x = =bxo, and tilted 
at the angle of it arctan(|c|). 

We now consider the one-dimensional critical NLS (16a) with the two tilted-beams initial con- 
dition 

^o(^) = (1 - ^)Citcit(0' x-xo) + {l-{e + Ae))^*;'*-[,(0, X + xo))e'^', (42) 

where V'o ei^^ defined in (41) and e = 10"'^. This initial condition correspond to two input beams, 
centered at ixo, titled toward each other, possible with a different power (when Ae ^ 0), and with 
a relative phase difference A9. In Figure 4(a) we plot the solution when Ae = and A9 = 
(equal-power, in-phase input beams). Since the power of each beam is slightly below Per, each 
beam focuses up to a certain time {t w 0.3), and then defocuses. Subsequently, the two beams 
intersect around t ~ 0.6. Since the beams are in phase, they interact constructively. As a result, 
their total power is w 2Pcr- Hence, the solution collapses at Tc « 0.6. 

In Figure 4(b) we repeat this simulation with Ae = and AO = vr (equal-power, out-of-phase 
input beams). Before the two beams intersect, their dynamics is the same as in Figure 4(a). When 
they intersect at t ~ 0.6, however, the two beams are out-of-phase. Hence, they repel each other. 
Since each beam has power below Per, there is no collapse. 

In Figure 4(c) we repeat this simulation with AO = and Ae w 2.7- 10~^ (in-phase and slightly- 
different input powers), and observe that at i ~ 0.6 the two beams repel each other and there is no 
collapse. In particular, comparison of Figures 4(a) and 4(c) shows that the 0(10"^) change in the 
initial condition lead to a completely different "post collapse" interaction pattern between the two 
beams. 

The dynamics in Figure 4(c) is qualitatively the same as in Figure 4(b). This suggests that 
when the two beams in Figure 4(c) intersect, their phase difference is ~ vr. Indeed, let il:i{t,x) 
and V'2(i) x) be the solutions of the one-dimensional critical NLS (16a) with the initial conditions (1— 
^)V'expikiit(0'^ ~ ^o) and {I - {e + Ae))^**J'iri^.(0, x + xq)), respectively. In Figure 5 we plot the 
difference between the phases of ■0i(i, x) and 'i/'2(i,3;), and observe that around t 0.6, this phase 
difference is indeed ~ vr. 

We thus see that. 

Conclusion 1. Because of the loss of phase after the collapse, the phase difference between post- 
collapse intersecting beams becomes unpredictable. 

Therefore, as noted by Merle [2], the interactions between two post-collapse beams are chaotic. 
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Figure 4: Solution of the one-dimensional critical NLS (16a) with the initial condition (42) with Tc = 0.25, xo ~ 6, c = 
8 and e — 10^'^. (a): In-phase identical beams {A9 — and Ae — 0). (b): Out-of-phase identical beams {A6 = tt 
and Ae — 0). (c): In-phase non-identical beams {A9 = and Ae ~ 2.7 ■ 10"'*). 
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Figure 5: a.Tgtpi{t,x = xq — ct) — a.rg^2{t,x = —xq -\- ct). 



8. Ring-type singular solutions 

In this section we propose a continuation of ring-type singular solutions, which is based on 
adding a reflecting hole with radius tq around the origin, and then letting tq — t- 0+. 

8.1. Theory Review 

Consider the two-dimensional, radially-symmetric critical NLS 



r 

Let us denote the location of the maximal amplitude by 

rraax{t) = argmaxjV'l. 



(43) 



(44) 
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Singular solutions of (43) are called 'peak-type' when rmax{t) = for < t < Tc, and 'ring-type' 
when TjYiax 

(t) > for < t < Tc. 

Let 



where 



L{t) = y/l-aH, T = [ Y^ds = ln(l - a^t), (45b) 

Jo ^ (s) " 

and G{p) is a solution of 



r" 

G"(p) + - + 
P 



4 

16^ 



G + g3 = 0, 0/G(0)gR, G'(0)=0. (46) 



In [24], Fibich, Gavish and Wang showed that i/jq is an explicit ring-type solution of the radially- 
symmetric critical NLS (43) that blows up (in L^) at Tc = Setting t = in (45) gives the 
corresponding initial condition 

V;^(r) =G(r)e^*(°'/^)'^'. (47) 

Equation (46) has the two free parameters a and G(0). However, in the case of a single-ring 
G profile, these two parameters are related [25]. For example, in the numerical simulations in this 
section, the G profile is the single-ring solution of (46) with 

G(0) PS 7.6 • 10"^ a ^ 0.357832. (48) 

8.2. Vanishing-hole continuation 

The sub-threshold power continuation approach of Section 4 cannot be applied to ipQ^\ since 
these solutions have an infinite power. In addition, this continuation cannot be applied to the 
ring- type singular solutions which are in H^, since these solutions exist only for P ^ Per [24]. 
Therefore, we now develop a different continuation approach, which is based on a vanishing-hole 
limit. 

Let 'ip{t,r) be a shrinking-ring singular solution of the critical NLS (43) with an initial condi- 
tion ifjoir). Let us add a hole around the origin with radius tq, and impose a Dirichlet boundary 
condition at r = rg, which is equivalent to placing a reflecting conductor at r = vq. In order for 
the initial condition ^^(r) to satisfy the Dirichlet boundary condition, we slightly modify it with 
a cut-off function Hs{r/rQ), i.e.,^ 



V'&(r)^V^^''(r) ■.= i^Ur)-Hs 



We thus solve the two-dimensional, radially-symmetric critical NLS 

I 

iil^tit, r) + iprr + -V'r + j'^I^V' = 0, tq < r < oo, (49a) 
r 



For example, in our simulations we used the cut-off function 

< p < 1, 

8 V 2/ 2 




Hs{p)={ Psip) 1<P<2, , P,{p) = 6(p-^] -5(p-^] +^fp-^^ ^ 



P> 2, 
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with the initial condition 



r 



and the Dirichlet boundary condition 

ip{t,ro) = 0, t>0. (49c) 

A typical simulation is shown in Figure 6. Initially, the ring solution shrinks and becomes higher 
and narrower as it approaches the hole. After the solution is reflected outward by the hole, it 
expands and becomes lower and wider. 

We now show that the conditions of Continuation Result 2 hold: 

1. The solution of (49) exists globally, since otherwise it collapses at some < < oo (a 
standing ring), or at r — )■ oo (an expanding ring). The flrst possibility is only possible, 
however, if the nonlinearity is quintic or higher, and the second possibility is not possible for 
any power-nonlinear ity, see [26]. 

2. By continuity, limrp_>o+ ip = V'g^^ fo^^ < t < Tc. 

3. The solution of (49) is invariant under the transformation (38). 

4. limt^Tc arg^^^^ = oo, see (45b). 

5. Let rj["°r := ar g mint r max (t) denote the reflection time, where rmax{t) is given by (44). Since 



the solution focuses for < t < T^°f and defocuses for T^°r < t < oo, it is collimated 



" ref 

at t = T^^j.. Therefore, 

arg V(T;°^ ,r;ro) = , a'^ E R. (50) 
Hence, by the arguments in the proof of Continuation Result 2,^ 

^(r;»^ + t, r; ro) = e*^(^'o)r(C^ " t, r; ro), t> 0. (51) 

Indeed, in Figure 7 we plot the solution at t = Tj^r it At for three different values of At, and 
observe that in all three cases, the two curves lie on top of each other. 

Therefore, we get the following result: 

Continuation Result 3. Let iJj{t,r;ro) be the solution of the NLS (49), and assume that (50) 
holds. Then, for any G M, there exists a sequence ro,n — ^ 0+ ( depending on 6 ), such that 

hm Mtr-r. )-[ < t < T,, 

^^hm^V'(t,r,ro,.) - j ^*^(e.)^2T,-t,r)e^^ %<t, ^^'^ 

where 'il)Q^\t,r) is given by (45), and Tc = Xjc?. Hence, this continuation satisfies Properties 1 
and 2. 

In particular, the limiting width is given by 



lim L(t;ro) 



1-1 
Tc 



*Here T^°f is the analog of T^^^, and (51) is the analog of (39). 
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Figure 6: Solution of the NLS (49) with n, = 0.08 at (a): to « 7.806, (b): ti « 7.808, (c): t2 « 7.8095 , (d): ts « 7.8100, 
(e): ti ~ 7.8110, (f): ts « 7.8126. Tlie arrows denote the direction in which the ring moves. 



< = r;", ±0.003. 



t = r"v ±0.0039 
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Figure 7: Solution of Figure 6 at t = T^^j. - At (solid) and at t = T^"^ + At (dashes), where T^°j: « 7.8096. The 
dotted line is the best fitting profile. All three curves are indistinguishable. (a): At = 0.0035. (b): At = 0.0037. 
(c): At = 0.0039. 



8.2.1. Simulations 

In Figure 8 we plot the reflection time as a function of rg, and observe that this data is in 
excellent fit with the parabola 

Tref = Tc + k2 ■ rl fc = 7.80981, A^a = -0.03585. (53) 

I 'jp ^ I 

Since Tc = ^ 7.80983, see (48), the extrapolation error is ' "j, = 0.0002%. The observation 
that Tc — Tref{ro) scales as Tq and not only as tq, will allow us to extrapolate ipQ {t, r) and L(t; rp) 
as a function of Tq, rather than of ro, leading to more accurate extrapolations. 
We now present simulation results that support Continuation Result 3: 

1. In Figure 9(a) we plot L(i;ro) for ro = 0.15,0.125,0.1 and 0.08. The curve which is obtained 
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from the extrapolation of these curves to rl = 0, is nearly identical to the limiting curve L 



1 - ^ 



, see Figure 9(b). 




0.05 , 0.1 0.15 



Figure 8: The reflection time T^°f as a function of the hole radius ro, for the solution of the NLS (49). Solid line is 
the parabola (53). 
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Figure 9: Solution of the NLS (49). (a): Solution width for ro = 0.15 (dots), ro = 0.125 (dashes-dots), ro = 0.1 
(dashes) and ro = 0.08 (solid), (b): Extrapolation of the curves {L{t;ro)} from (a) to ro = (dashes). Solid line 

is L{t) = ^\i-t/n\. 

2. In Figure 10(a) we fix the time at ti = 7.8124 > Tc, and plot the solution profile for tq = 
0.15,0.125,0.1 and 0.08. The curve which is obtained from the extrapolation of the pro- 
files {\'ip{ti,r;rQ)\} to rg = is nearly identical to \'ipQ^^\2Tc — ti,r)\, see Figure 10(b). 
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Figure 10: (a): Solution of the NLS (49) at t = 7.8124 = l.OOOSTc for ro = 0.15 (dots), ro = 0.125 (dashes- 
dots), ro =0.1 (dashes) and ro = 0.08 (solid). Solid bold line is the extrapolation of these curves to ro — 0. (b): The 
extrapolated proflle from (a) (solid bold). Dotted line is |7/;Q'^^'(2rc — 7.8214, r)|. The two curves are indistinguishable. 
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3. In Figure 11 we plot the accumulated phase at the ring peak, i.e., arg ip'^^{t,rmax{t)), and 
observe that small changes in rg hardly affect the phase before the singularity, but lead to 0(1) 
changes in the phase after the singularity. 

120 
100 

arg v|/(t,0) 

80 
60 

7.8 7.805 I 7.815 

Figure 11: Accumulated phase as a function of t, for the solution of the NLS (49) with ro — 0.08 (solid), ro — 0.1 
(dashes), ro = 0.125 (dots), and ro = 0.15 (dashes-dots). 



9. Vanishing nonlinear-damping solutions 

In this section, we propose a continuation which is based on the addition of nonlinear damp- 
ing. The motivation for this approach comes from the vanishing-viscosity solutions of hyperbolic 
conservation laws. Of course, the key question is which physical mechanism should play the role 
of "viscosity" in the NLS. In the nonlinear optics context, there are numerous candidates, which 
correspond to the mechanisms that are neglected in the derivation of the NLS from Maxwell's equa- 
tions: Nonparaxial effects, high-order nonlinearities, dispersion, plasma effects, Raman, damping, 
etc. Of course, for a physical mechanism to be able to play the role of "viscosity", it should arrest 
the collapse regardless of how small it is (so that we can take the limit of this term to zero, and still 
have global solutions). This requirement rules out some candidates (such as linear damping, see 
below), but still leaves plenty of potential candidates (such as nonlinear saturation, see Section 6). 

In this study we consider the case when the role of viscosity is played by nonlinear damping. 
The addition of small nonlinear damping is "physical". Indeed, in nonlinear optics, experiments 
suggest that arrest of collapse is usually related to plasma formation, and nonlinear damping can 
be used as a phenomenological model for the multi-photon absorption by the plasma. In BEG, 
a quintic nonlinear damping term corresponds to losses from the condensate due to three-body 
inelastic recombinations. In [27], Bao, Jaksch and Markowich showed numerically that the arrest 
of collapse by the addition of a quintic nonlinear damping to the cubic three-dimensional NLS is 
in good agreement with experimental measurements. Nonlinear damping arises also in the context 
of the complex Ginzburg-Landau equation, see Section 10. 

9.1. Effect of linear and nonlinear damping - review 

In [28], Fibich studied asymptotically and numerically the effect of damping on blowup in the 
critical NLS, and showed that when the damping is linear, i.e., 

i?/'t(i,x) + A?/; + |V'|^/'^V' + ^<5V' = 0, V(0,x) = V'o(x), (54) 

if the initial condition ipo{x) is such that the solution of (54) becomes singular for 6 = 0, then the 
solution of (54) exists globally only if 5 is above a threshold value 6c > (which depends on i^q). 
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Therefore, linear damping cannot play the role of viscosity in defining weak solutions of the NLS. 
When, however, the damping exponent is critical or supercritical, i.e., 

iV'i(t,x) + AV + (l + i<5)|V|^/'^V = 0, <5>0, (55) 



or 



z^j(t,x) + + IV'l^/'^V + ^'^lV'rV' = 0, (5>0, p>4/d, (56) 



respectively, then regardless of how small 5 is, collapse is always arrested. Therefore, Fibich sug- 
gested that nonlinear damping can "play the role of viscosity" in defining weak NLS solutions, i.e., 
we can define the continuation 

lb := lim 'ib^^\ (57) 
5^0+ 

where ip^^^ is the solution of (55) or (56). 

Since the results in [28] are not rigorous, we now present the relevant rigorous results that exist 
in the literature. Passot, Sulem and Sulem proved that high-order nonlinear damping always 
prevents collapse for d = 2. Antonelli and Sparber extended this result to d = 1 and d = 3: 

Theorem 7. [29, 30] The d- dimensional cubic NLS with nonlinear damping 

zV't(t,x) + + A|^/'|V + «'5|';/'r~V = 0, AgM, (5>0, (58) 

where ^o(x) G H'^iW^), 3 < p < oo if d = 1,2, and 3 < p < 5 if d = 3, has a unique global in-time 
solution. 

This rigorously shows that high-order nonlinear damping can play the role of "viscosity". More 
recently, Antonelli and Sparber proved global existence for the case where the damping exponent 
is equal to that of the nonlinearity: 

Theorem 8. [30] Consider the cubic nonlinear NLS with a cubic nonlinear damping 

iV'i(t,x) + AV' + (l + i<5)|^|V = 0, (59) 

where V'o(x) G H^{W^), x^/^o ^ L'^{W^), and d < 3. Then, for any 5 > 1, equation (59) has a unique 
global in-time solution. 

Unfortunately, because of the constraint 5 > 1, Theorem 8 does not show that critical nonlin- 
ear damping can play the role of viscosity. We note, however, that the asymptotic analysis and 
simulations of [28] strongly suggest that the solution of (55) exists globally for any < 5 <C 1. 

9.2. Explicit continuation when ipoir) = ^explicit (0, ?") 

In the special case where tp^^^ is the solution of (55) with the initial condition ^explicit (^ = 0), 
we can calculate explicitly the vanishing nonlinear-damping limit (57): 

Continuation Result 4. Let ip^^\t,r) be the solution of the NLS (55) with the initial condition 

ipo (r) = '(/'explicit (0, r ) , (60) 
see (7) . Then, for any G M, there exists a sequence 5n — >• 0-1- ( depending on 6 ), such that 

^ - 1 V':xp,icit,.(2T, - t, r)e^^ < t < oo, ^^^^ 
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where ^cxpiicit.K is given by (8a) with a = k, 

K = 7r [A(0)4(^*) - MO)Bi{s*)] « 1.614, (62) 

Ai{s) and Bi[s) are the Airy and B airy functions, respectively, and s* ~ —2.6663 is the first negative 
root ofG{s) = V2>Ai{s) - Bi{s). 

In particular, the limiting width of the solution is given by 

lim L{t; 5) = | \~ % , IV (63) 
5^0+ ^ ' [ K{t-Tc) Tc<t <oo, 

Proof. See section 9.4. □ 

Equation (61) provides a continuation of the explicit blowup solution '(/'explicit beyond the sin- 
gularity. As with all the continuations that we saw so far, the weak solution is only determined 
up to a multiplicative phase constant e*^ (Property 1). In contrast with these continuations, how- 
ever, the continuation (61) is asymmetric with respect to Tc, since k ^ 1. Therefore, it does not 
satisfy Property 2. This is due to the directionality in t of the damping effect, i.e., the fact that 
equation (55) is not invariant under the transformation (35). 

Remark 1. The value of k is also given by, see (94), 

A{0) 



-Ai{s*) 

Remark 2. The limiting solution in Continuation Result 4 does not have a non- collapsing "tail", 
since the power «/ V'cxpiicit k(^-^c — t,r)e^^ is equal to that of 'iJjo{r). 

Remark 3. In order to understand why Li increases (and not decreases) after the singularity, we 
note that while the vanishing nonlinear- damping does not affect the solution power, it increases the 
Hamiltonian, see Section 9.5. Since 

I? 

i?(V'expiicit,a) = -jWR^'^W = Ma^ 

the increase in the Hamiltonian implies that the defocusing velocity (angle) should be higher than 
the focusing velocity (angle), see Figure 12. 

9.2.1. Simulations 

In order to provide numerical support to Continuation Result 4, we solve numerically the damped 
NLS (55) with d = 1 and the initial condition (60) with Tc = 1, for various values of 5. Figure 13(a) 
shows that as 5 — )■ 0-|-, the maximal amplitude increases, and that it is attains at Z^a^. — ?• T^.. 
Figures 13(b) and 13(c) show that lim5_^.o+ L{t; 5) is given by (63), and that 

5^0 I ^ Tc<t <co, 

respectively, where k is defined in (62). 

In Figure 13(d) we plot the accumulated phase at x = 0, and observe that small changes in 5 
have a negligible effect on the phase before the singularity, but an 0(1) effect on the phase after the 
singularity, which is an indication that phase of the weak solution becomes non-unique for t > T^- 
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Figure 12: The vanishing nonhnear damping hmit (61) for d = 1 and Tc = 1. (a) Color plot, (b) Contour plot, (c) 
Surface plot. 
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Figure 13: Solution of the damped NLS (55) with d = 1 and the initial condition (60) with Tc = 1 for 5 = 10~^ 
(dashes-dots), & = 10~^ (dashes) and 5 = 1.25 • 10~^ (dots), (a): ||V'l|oo. (b): L{t;5), recovered from ip using (17). 
Solid line is (63). (c): Lt{t;5). Solid line is (64). (d): Accumulated phase at a; = 0. 



9. 3. Continuation for loglog collapse 

Let be a solution of the undamped critical NLS that undergoes a loglog collapse. Since 
nonlinear damping leads to defocusing (and not to oscillations) after it arrests the collapse [28] , the 



24 



limiting solution has a point singularity and not a filament singularity. Indeed, the continuation 
has an infinite-velocity expanding core:^ 

Continuation Result 5. Let ipo{r) be a radial initial condition, such that the corresponding so- 
lution of the undamped critical NLS (5) collapses with the iI'rW profile at the loglog law blowup 
rate at Tc- Let il)^^^ be the solution of the damped NLS (55) with the same initial condition. Then, 

lim V^'^^ = V^, < t < Tc. 
(5-5.0+ 

In addition, for any < 5 ^ 1, there exists 9 (6) G M, and a function (p £ L'^ , such that 



lim 

<5^0+ L 



t n+, 



where ipjim is given by (25b) with some function L{t;5), such that 



lim lim L(t:5)=0, lim lim Lf(t: 5) =oo, lim 9(6) = oo. 

Proof. See Section 9.6. □ 

The post-collapse infinite velocity of the expanding core, is a consequence of the infinite velocity 
of the loglog collapse before the singularity, and the increase of the velocity after the singularity 
(see Remark 3). 

Remark 4. Because of the infinite velocity of the expanding core, it "immediately" interacts with 
the non- collapsing tail. Therefore, the validity of the reduced equations that are used in the deriva- 
tion of Continuation Result 5 breaks down "shortly" after the arrest of collapse. See Section 9.6.4 
for further discussion. 

9.3. L Simulations 

In order to illustrate Continuation Result 5 numerically, in Figures 14-15 we solve the damped 
NLS (55) with = \/r05 V'expiicit (* = 0) for 5 = 2 • 10-^ 6 = 1.5 ■ 10'^, and 6 = W^, and 
observe that the solutions are highly asymmetric with respect to Tmax- In addition, as 5 — )■ 0+, 
the post-collapse expansion of the singular core becomes faster and faster. 

9.4- Proof of Continuation Result 4 

In order to prove Continuation Result 4, we first approximate the NLS (55) with a reduced 
system of ordinary differential equations. Then, we solve the reduced system explicitly as (5 — t- 0+. 

9.4-1- Reduced equations 

Lemma 3. Let ip^^^ be the solution of the damped NLS (55) with the initial condition (60). Ifip^^^^ 
i^RW) see (25), the reduced equations for L[t) are given by 

m = -^, ~S= (65a) 
Luit) = -^, (65b) 



^The observation that the velocity of the expanding solution is infinite, is due to Merle [31] 
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Figure 14: L(t) for the solution of the damped NLS (55) with d = 1, and the initial condition '0o(r) 
y055V)expiicit(t = 0,r;Tc = 1) for & = IQ-^ (solid), S = 2- 10"^ (dots), and 6 = 2.5- 10"^ (dashes). 
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Figure 15: Same as Figure 12 for the numerical solution of the damped NLS (55) with d = 1, 5 = 10 ^, and the 
initial condition ipo{r) = VTOS V'expiicit (i = 0,r;Tc = 1). 



where q = 1 1 -R^^^ 1 1 4/^+3' ^ 

is given by (27), with the initial conditions 



m = 0, L(0) = Te, Lj(0) = -1. 



(66) 



Proof. In [28], Fibich used modulation theory [9] to show that when V'corc ~ '4'r(o), self focusing 
in (55) is given, to leading-order, by 



Pt{t) - — ^2 ivri2' ~ ~ 



(67) 



where v{(3) is given by (27). We recall that when 5 = 0, the initial condition (60) leads to the 
explicit blowup solution (7), for which L{t) = Tc — t. Therefore, when 6 = 0, 



L(0) = T„ Lt{0) = -1, Lu{0) = 0. 
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The initial condition is independent of the subsequent dynamics, hence it is independent of 5. 
Therefore, the initial condition is also given by (68) for 6 > 0. Therefore, since /3(t) = —L^Lu, 
then /3(0) = 0. Now, since w(/3 = 0) = 0, see (27), by (67) we have that pt < 0. Therefore, /3 < 0, 
and consequently v{l3) = 0, see (27). □ 

9.4-2. Simulations 

The derivation of the reduced equations (65) is based on modulation theory, which is not rigorous. 
Therefore, we now provide numerical support for the validity of (65). In Figure 16 we solve (65) 
numerically for d = 1 and various values of 5. We compare these solutions with direct simulations 
of the damped NLS (55) with d = 1 and the initial condition (60), from which we extract the value 
of L{t) using (17). When 6 = 10"^, the two curves are similar, and for S < W ^ the two curves 
are indistinguishable. This confirms that as (5 — )• 0+, the dynamics of the solution of the damped 
NLS (55) with the initial condition (60) is given by the reduced equations (65). 

In Figure 17 we plot the rescaled profile L^/'^{t)\ij{t,x/L)\ at t = 0.6 < and at t = 1.4> T^. 
The two rescaled profiles are indistinguishable from each other and also from R{x), providing 
support to the assumption that ~ i^RW , which was used in the derivative of the reduced 
equations [9]. 




0.5 1 1.5 0.5 1 1.5 

t t 



Figure 16: Width L{t) of the solution of the damped NLS (55) with d = 1 and the initial condition (60) with Tc — I 
(solid). Dashed line is the solution of the reduced equations (65)-(66). (a): S = 10~^. (b): S = 10"'*. (c): S = 10~^. 
(d): S = 1.25 • 10"^ 



9.4.3. Analysis of the reduced equations 

Our ultimate goal is to solve the ODE system (65) with the initial conditions (66) explicitly 
as 5 — )• 0+. We first prove the following Lemma: 
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Figure 17: Rescaled profile L^^'^{t)\il}{t,x/L)\ of the solution of the damped NLS (55) with d = 1 and the initial 
condition (60) with Tc = 1 and 5 = 1.25 • 10"'^, at t = 0.6 < Tc (solid) and at t = lA> Tc (dots), where L{t) is given 
by (17). The dashed curve is R{x). All three curves are indistinguishable. 



Lemma 4. The solution of the reduced equations (65) with the initial conditions (66) can he written 



as 



L{t) = 1/Aisit)), 



where 



A{s) = TT 



V3Ms)-Bi{s) 



+ 



1 Bjjs) 



Ai{s) and Bi{s) are the Airy and B airy functions, respectively, and 

s{t) = -5^/3 f A^{w)dw. 
Jo 

Proof. Equation (65a) can be rewritten as 

Since = 0) = /3(r = 0) = 0, 
We recall that [9] 



(69) 



(70) 



1 



L^is] 



ds. 



Att 



-6t. 



A 



1 



(71) 



(72) 



Substituting (71) in (72) gives Ad- 
equation into Airy's equation 



A ' L 
-6tA. The variable change s = —5^^^t transforms this 



A.s.s — sA. 



(73) 



Since A = 1/L, then by (66) and (70), the initial condition for (73) are 



A{s = 0) = A{t = 0) 



T' 



(74a) 

As{s = 0) = -5~^/^At{t = 0)A-'^{t = 0) = -6~^/^. (74b) 
The solution of Airy's equation is a linear combination of the Airy and Bairy functions: 

A{s) = hAis) + k2Bi{s). (75) 
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Substituting (74) in (75) gives 



~§-i/SB.io) + ±Biio) _ rV3A,(o) + ^a:(o) 



= . rrhr^rTTTT, ko 



The Airy and Bairy functions satisfy the Wronskian relation 

Ais)B',{s)-A',is)B,is) = -. (76) 

TT 

Therefore, 

fei = vr (^6-'/'BM + ^^KO)) , k2 = -n (^6~'/''M0) + 1a^(0)) . (77) 

The Airy and Bairy function satisfy 

MO) = B,{0)/V3, MO) = -B'M/VS. (78) 
Substituting of (77) in (75), using relations (78) and (76) leads to (69). □ 



9.4.4. S-^0+ 

Lemma 5. Let be the first negative root of A{s;5), see (69). Then, \Ynis^s*_ t(s;6) = 00. 
Proof. Inverting (70) gives 





Let < e < 1. Then, 



By definition, A{s'~;6) = 0. Furthermore, since A{s;6) is a nontrivial solution of Airy's equation, 

then As(s|;(5) 7^ 0, since otherwise from uniqueness it follows that A(s;5) = 0. Therefore, there 
exists < e such that 

A{s;~5)^{s-s*g)As{sf,~5), s| < s < 5? + (79) 

Hence, 



A^s) Ajisl,6)Jst {s-s*-)^ 



Therefore, t{.s*^) = 00. □ 



By (70), s{t) is monotonically decreasing from s(t = 0) = 0. Hence, the interval < t < 00 
transforms to > s > s^. Since the Airy and the Bairy functions are bounded for s < 0, then A{s; 6) 
is finite for s < 0, see (69). This shows that the solution of the damped NLS (55) does not collapse. 
Note that A{s = sp = corresponds to an infinite beam width, i.e., to a complete diffraction. 
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Lemma 6. Let s* = lini^_j_Q_|_ si. Then, s* ~ —2.6663 is the first negative root of 

G{s) = 0, G(s) := VsAiis) - Bi{s). (80) 
Proof. By (69), as 5^ 0+, 

A{s) ~ 7rAi(0)<5~i/3 \V3Ai{s) - Bi{s)\ = 7r^i(0)r ^/^^(s). (81) 

Therefore, s* satisfies (80). The value of s* was computed numerically, see Figure 18(a). □ 




(c) 



1.5 -3 s* -2 „ -1 



Figure 18: (a): The function G(s), see (80). Here, s* is the first negative root of G{s). (b): s{t), as calculated from 
numerical integration of (70) for S = 10"" (dashes-dots), 5 = 1Q~'^ (dashes), 5 = 10"*^ (dots) and 5 = 1.25 • lO"'' 
(solid), (c): Same as (b) for the inverse function t{s). 

In Figure 18(b) we plot s{t) by numerically integrating (70), where A(s) is given by (69), and 
observe that as 5 — )■ 0+, s{t) tends to the step function 



lim s{t) 



0<t<Tc, 
Tc<t<oo. 



Therefore, the inverse function t{s) tends to the step function. 



lim t{s) 



oo s = s 







s* < s < 0, 
s = 0, 



see Figure 18(c). We now prove these limits analytically. 

Lemma 7. For any si such that —1 ^ si < 0, lim^_j,Q_|_ t(si) = Tc. 

Proof. Equation (69) can be rewritten as 

A{s;~5) = R{sy5~'/'' + ^^Fis), 



where 



R{s) = nAi{0)G{s), F{s) 



Bijs) 



TrA'MGis). 



(82) 
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-ds 



^1/3 



R{s)+'^F{s 



R{0) = 0, F(0) = 1, 
R'{0) = irAiiO) (^3^(0) - ^,'(0)) = -27rAi(0)S,'(0). 



Therefore, 

t{si) = 5-1/3 

By (78), (80) and (82), 
and 

Furthermore, by (78), 

MO)Bi{0)+A'i{0)Bi{0) = 0, 
and by the Wronskian relation (76), 

Adding (86a) to (86b), gives 27rAi(0)5,'(0) = 1. Therefore, by (85), 

R'{0) = -1. 

Thus, by (84) and (87), 

R{s) s, F(s) ~ 1, -Ks<0. 

Substituting this in (83) gives 



tcIs. 



Therefore, hm^_^Q^ i(si) = Tc. □ 
From (88) it follows that 



51/3 



-s + 



51/3 



51/3 



51/3 



-si+ 



lim t{s = -c^i/3) = T, 



(83) 

(84) 
(85) 
(86a) 

(86b) 
(87) 



(88) 



Therefore, t{s) has a boundary layer at s = with thickness 5^^^, in which it increases monotonically 
from to Tc- Hence, 

Corollary 4. lim^^g_^ s{t; 6) = 0, 0<t<Tc. 
Lemma 8. For any s* < S2 < si < 0, lim^^Q {t{s2) — t{si)) = 0. 
Proof. By (83), 



tis2)-t{si) 



51/3 



S2 



R{s) + ~^F{s) 



Since R{s) = 7rAj(0)G(s) and G{s) > in (s*, 0), then R{s) > c> for any s* < S2 < s < si < 0, 
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where c is independent of 5. Therefore, as 5 — )• 0+, 



hm(t(s2) -t(si)) ~ hm(5^/3 / 

5^0 <5-s>0 Js2 



ds = 0. 



□ 



From Lemma 7 and Lemma 8 it follows that 
Corollary 5. 



lim t{s) = Tc, s* <s < 0. 
5-s>0 



Lemma 9. lim^_^Q^ s{t; 6) = s* for Tc < t < oo. 



Proof. Since s(t;6) is monotonically decreasing (see equation (70)), t{s;6) w Tc for s| < s < 
(Corollary 5), and lim5_^c,i t{s; 5) = oo (Lemma 5), then near s| there is a boundary layer in which t 
changes from to oo. Therefore, the values (T^, oo) are attained in the boundary layer around s|. 
Since lim? , s* = s*, the result follows. □ 

From Corollary 4 and Lemma 9 it follows that, see Figure 18(b), 
Corollary 6. 

hm sit-d) = \ T~<t< ^ ^^^^ 

(5->o+ y ^ ic < I < oo. 

Our next observation concerns L{t). 

Lemma 10. 

5^0+ [ K{t-Tc) Tc < t < OO, 

where k is defined in equation (62). 
Proof. Using L = A^^ and (75), 

L, = -^-'^^ = -A-'A,{-~5'I^A') = ~5'I^A, = ~5'/' [k.A'M + hBUs)] . (91) 

By (77), 

lim 5^/3^1 = TTBi{0), lim ^^/^fes = -7rAi(0). 
5^0+ 5^0+ 

Therefore, 

lim Lt{t) = TT [B,{0)A'Mt)) - Ai{0)Bi{s{t))] , s{t) = lim s(t; ,5). 

5^0 S^O 



Hence, by (89), 



limLt(t) 

5^0 



TT [Bi{0)A'.{0) - Ai{0)B',{0)] < t < Tc, 
TT [Bi{0)A[{s*) - A,{0)Bi{s*)] Tc < t < oo. 



Since vr [Bi{0)A'i{0) - Ai(0)5|(0)] = -1, see equation (76), 



limLt(t) 

5^0 



-1 < t < Tc, 
Tc < t < oo. 
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where k is given by (62). Since L(0) = Tc, (90) follows. □ 



Our last observation concerns the solution phase. By definition (25b), 



* 1 



Therefore, by (90), 



;im T(r,) = / ' ^ ds = oo. (92) 

5->0+ Jo (?c - SY 



Therefore, the phase information is lost at the singularity. 

9.4-5. Proof of Continuation Result 4 

We have that i/j^ ~ V'_r(o); where L(t) is given by (90). Therefore, by Lemma 2, when < t < 
Tc, '4'R{o){r,t) = V'explicit(i,'r), and when Tc < t, ipji(o)(t,r) = 'i/'explicit,a(2^c - t,r). In addition, 
since arg'0(t,O) ~ arg ^^(o) (t, 0) = r(t), equation (92) shows that the limiting phase becomes 
infinite at and after the singularity. Hence, for a given given t > Tc and G M, there exists a 
sequence 6n 0+, such that lim5^_^o+ arg V'^^^H^i 0) = ^- Hence, Continuation Result 4 follows. 

9.5. Hamiltonian dynamics 

In the case of a non-conservative perturbation such as nonlinear damping, the Hamiltonian of 
can be approximated with, see [9, eq. (H.5)], 

H{^) ~ ^{L\,. (93) 



Since = -^L\L%tt, then 



Therefore, the Hamiltonian increases with t. Moreover, by (70) and (81), 

Hs = = m4j^, = -Mp/'A^ ~ -M7r^AU0)G\s). 



Therefore, 



Now, 



AH := H{t = oo) - H{0) = Hsds = M^^^^^O) / G'^{s) ds. 

Js=0 Js=s* 











G^{s)ds = G^s -/ 2sGGs{s)ds = -2 j sGGs{s)ds. 
Since G„ = sG, then 







s=s' 



Therefore, 



r-O 

G2(s) ds = -2 I GssGsis) ds = -Gl 



AH = -M7r2A2(0)G2 ° 



33 



Since the Wronskian of the Airy equation is a constant, 



W{G,A) = Gis)A'^{s) - Ms)G'{s) = GiO)A'M - Ai(0)G'(0). 
Therefore, since G vanishes at s = 0,s* , 

Also, by (87), 

-1 = R'{0) = 7rAi(0)G'(0). 

Therefore, 

2 



AH = -MTr^Af{0)Gi{0) 1 



Ai{0) 




Ms*). 

On the other hand, by (93), H ~ M(Lf - P/L'^). In addition, since the limiting solution has 
exactly the critical power, then lim^_^Q^ /3 = 0, see (28). Therefore, lim.^_^Q^H = ML^. Hence, 
by (64), 

AH = M{k^ - 1). 

Comparison of the last two expressions for AH shows that 

AM 



\Ai{s* 



(94) 



9. 6. Proof of Continuation Result 5 
9.6.1. Analysis of the reduced equations 

We first analyze the dynamics of the collapsing core under the assumption that it is governed 
by the reduced equations (67).^ By continuity, as 5 — > 0+, the limiting solution undergoes a 
loglog collapse as t — > Tc—. Therefore, as t — > Tc—, the amount of power that collapses into the 
singularity is exactly Per- Hence, by (28), lim5_j.o+ P{Tc) = 0. Therefore, since f3t < 0, see (67), and 
since i^(/3) = for /3 < 0, we have that lim5_5.o+ /3f = for t > Tc. Therefore, lim5_5.o+ /3 = 0— for 
t > Tc- Hence, the collapsing core expands linearly with the ip^iio) profile. Therefore, the expanding 
core is given by ^cxpiicit,a(2'^c -t,r). 

The above arguments show that according to the reduced equations, after the singularity the 
solution is of a Bourgain-Wang type, but do not provide the value of the expansion velocity a. In 
order to do so, we now solve the reduced equations. By (67), 

/3r = -6, T = 1^ 12^"^'- 

For solutions that undergoes a loglog collapse we have that /3(0) > 0, see (28). For a fixed 6 > 0, 
as /3 \ 0+, i/(/3) becomes negligible compared with 6. Therefore, to leading order, /3r = —5."^ 
Hence, 

/3(r) = /3o - 6t, /3o = /3(0) > 0. (96) 



® The validity of this assumption is discussed in Section 9.6.4. 

^ In other words, we approximate (67) with (65). The vahdity of neglecting v{l3) is discussed in Section 9.6.4. 
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Substituting (96) in (72) gives Arr = (/3o — 6t)A. The variable change 

s = r2/3/3o - S^y^T (97) 

transforms this equation into Airy's equation (73). 

Let A{t = 0) = ^0 = 1/^0 and At{t = 0) = A'q = -L'q/{LqY be the initial conditions for A{t). 
Therefore, the initial conditions for A[s; 5) are 



As{s = so) = -ri/3^t(t = o)A-2(t = o) = -rV3_o_^ (ggb) 



^(■5 = So) = -4o, (98a 
where 

so := 5-2/3 /3o. (98c) 

Therefore, A{s; 5) is given by (75) with 

fci = vr (^^r i/3^,(so) + AoBliso)^ , k2 = -7T (^^5-'/^A,iso) + ^oA-(so)) • (99) 

9.6.2. 6^0+ 

In Figure 19(a) we plot t{s;5), and observe that 

Lemma 11. 

C s = So, 

lim t (s; S) = { Tc s* < s and sq - s :^ 1 (100) 

-^^0+ ^ ^ oo s = s*. 



t{s;S) 




sit 6) 





(b) 







Figure 19: (a) t{s-S). (b) s{t-5). For < 5 < 1. 



Proof. By (95) and (97), ^ = -6^/^A^. Therefore, 

1 . 1 



t(s;5 



dt 
1 



51/3 4 A^s 



-ds 



51/3 [A;i^i(s) + A;25i(s)]' 



rds. 



(101) 



Following [32], it is easy to verify by direct differentiation and using the Wronskian relation (76) 
that 



t{s:5 



/fci(5V3 



Bijso) Bijs) 
A{so;6) A{s;6) 



IT 



Ai{sQ) Ai{s) 



A[so]5) A[s;5) 



(102) 
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Therefore, by (98b), 
t (s:6 



vr 



Ao A{s;6) 



vr 



Ai{s) A, (so) 



A{s;d) Ao 



(103) 



Lemma 12. Let s| be the first negative root of A{s;6). Then, Unisys*- iis',S) = oo. 
Proof. This foUows by direct substitution of s| in (103). □ 



Lemma 13. Let s* := hm^^j^Q^s?. Then, s* ~ —2.338 is the first (negative) root of Ai{s) 
Proof. By (98c), sq > 1 for < 5 < 1. Therefore, 



Ai{so) 



1 -1/4 _2s3/2 n / N 1 -1/4 2^3/2 
-Sn ^ e 3^0 , Bi{so) ^Sn ^ 63^0 . 



vr 



Hence, by (99), 



Therefore, 



lim — = 0. 



A{s;5) ~ A;i^i(s) 



5 = 0(1), 0<5<^1. 



Thus, the result fohows. □ 



From Lemma 12 and Lemma 13 it follows that: 
Corollary 7. t{s*;5) := lim^^^^ t(s|; ^) = oo. 
Equation (103) can be written as 



t s-5]=T 



1 



Ais;S) 
Aa J 



T, 



{5) _^ TT Bi{so) 

fciJVs Ao 



By (99) and (106), 



TTBi{so) 
fci .51/3^0 



(Ao)25i/3 



t:\t^^,+AoB[{so) 



(104) 



(105) 



(106) 



By (98c) and (104), B[{sq) ~ Bi{sQ) ■ s]^^ = Bi{so) ■ (J-i/^/?^^ Therefore, since L = 1/A, 



(Ao 



which is the adiabatic approximation of Tc, see [9, eq. (3.31)]. 

By (104), Bi{s) / Bi{so) is exponentially decreasing as s decreases from sq. Let us assume for 
simplicity that A'q > 0. Then, A{s;6) > A{so). Therefore, ^^(1'^^!^^^°^^ is exponentially decreasing 
as s decreases from sq. Hence, 

t (s; 6) ^ r]*^) w Tc, s* < s and So - s > 1. 
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This concludes the proof of Lemma 11. □ 



Prom Lemma 11 it follows that 

s{t) ~ I ^* rj^^^^ as 6^0+, (107) 

see Figure 19(b). We recall that 

Ltit; ~6) = 5^/3 [k.A'Mt)) + k2B'Mt))] , 

see (91). Hence, by (107), 

Lt{t;6)^S^/^[hA'i{s*) + k2Bi{s*)], T,<t, 0<5<1. 
Therefore, by (105), 

Lt{t; 6) ~ 6^/^hA'i{s*), n < t, < 5 < 1. 

Hence, L{t) is linear for < t. 
By (104) and (106), 

; H/3 ^ 1 n ^ \ 1 -1/4 2^3/2 

^0 -t c ^0 -t c 

Therefore, 



1^-1/4^.3/2^,,^, 

Since lim^_j^Q_|_ sq = OO; 
Therefore, the post-collapse slope of L[t) becomes infinite as 5 ^ 0+. □ 



V^^o e3^° T,<t, 0<(5«1. (108) 



lim 5^/3^i4(s*) = oo. 
(5->0+ 



9.6.3. Simulation 

In Figure 20 we solve the reduced equations (65), and observe that for t > Tc, Lt{t) is indeed in 
excellent agreement with the asymptotic prediction (108). 



91.37 
10 

L,(t) 



0.5 ^ 1 1.5 

Figure 20: Solution of the reduced equations (65) with d — 1, 5 = 2.5 ■ 10^*, and the initial conditions L{0) — 
1, L'q — —1, and /3(0) = 0.1. Sohd line is Ltit). Dashed line is the prediction (108). 
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9.6.4- Validity of the reduced equations? 

In the asymptotic analysis in Section 9.6 we approximated the critical NLS with the reduced 
equations (67). Then, we approximated the reduced equations (67) with the reduced equations (65), 
by neglecting i^(/3). We now consider the validity of these approximations. 

In Figure 21(a) we plot L{t) for the solutions of the reduced equations (67) and (65). Al- 
though 5 ~ 0.0019 is not much larger then z^(/3(0)) ~ 0.0016, the two solutions are "close". Plot- 
ting Lt{t) shows that in both cases, Lt{t) is a constant after the collapse is arrested, see Figure 21(b). 
The "addition" of i^(/3), however, decreases this constant by a factor of ~ 4. Therefore, when z^(/3) 
is not neglected, the approximation (108) for Lt{t) is not accurate, but the solution still expands 
linearly at a velocity that goes to infinity as 5 — t- 0+. 

X 10" 

61 ' ' ' 1 



L,(t) 
' 2 



: 

-2 1 ' ' ^ 

0.2 0.4 0.6 0.8 1 0.6 0.65 0.7 0.75 0.8 

t t 

Figure 21: Solution of the reduced equations (67) [solid] and (65) [dots] with S = 10^^, and the initial condi- 
tions ^(0) = 0.1, Lo = 1, and Lt(0) = -1. (a) L{t). (b) Lt{t). 

In order to check the validity of the reduced equations (67) with i^(/3), we solve the damped 
NLS (55) for d = 1 with the initial condition \/L05 V'expiicit (* = 0,r;Tc = 1), for various values 
of S. Then, we extract from ip the value of L(t) using (17). These NLS solutions are compared 
with solutions of the reduced equations (67). The initial conditions for the reduced equations 
are as follow. By (28), /5(0) ~ W^/^^e.pHatWl-P^r ^ q^o^ ^ 0.3242. By (17), L(0) ^ 0.9524. 
For ^cxpiicit(i = 0, r;Tc = 1), Lt{0) = —1, see (66). The multiplication by \/1.05 leads to a small 
change in Lt{0). We found that Lt{0) = —1.02 provides the best fitting. 

Figure 22 shows that there is a good agreement between L{t) of the reduced equations (67) and 
of the NLS. In addition, in both cases, the post-collapse defocusing velocity increases as 6 decreases. 
The curves of Lt{t) show a good agreement when the solutions focus, but differ when the solutions 
defocus. In particular, Lt{t) of the NLS solution is not a constant after the arrest of collapse. We 
relate this difference to the interaction between the expanding core and the tail, which is ignored 
in the reduced equations. This core-tail interaction did not occur in the explicit continuation case, 
see Section 9.2, since in that case the power of the initial condition i/'cxpiicit is equal to Per, hence 
there is no tail. In addition, this phenomenon did not occur in the sub-threshold continuation (see 
Section 4), since there the expansion velocity was finite. Therefore, sufficiently close to Tc+, the 
expanding core had a negligible interaction with the tail. 

In summary, the asymptotic analysis and numerical simulations suggest that in the nonlinear- 
damping continuation of NLS solutions that undergo a loglog collapse, the singular core of the 
NLS solution expands after the singularity with a velocity that goes to infinity as 5 — t- 0+. The 
post-collapse expansion velocity is, however, probably not linear in t. 
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0.45 0.5 0.55 0.6 
t 



Figure 22: Solution of the damped NLS (55) with d = 1 and the initial condition \/1.05'tp explicit {t — 0,r;Tc = 1) 
(solid). Dashed line is the solution of the reduced equations (67) with the initial condition /3(0) ~ 0.3242, -L(O) ~ 
0.9524, and Lt(0) = -1.02. (a) S = 10"^ (b) 5 = 1.5 • 10"^ (c) S = 2 ■ 10"^ Top row: L{t). Bottom row: Lt{t). 



10. Complex-Ginzburg Landau continuation 

The two-dimensional Complex Ginzburg-Landau equation (CGL) 

iil^tit, X, y) + A-0 + IV'PV' - ^eiV' - ie2'^ip + iealV'P^ = 0, 

arises in a variety of physical problems: Models of chemical turbulence, analysis of Poiseuille 
flow, Rayleigh-Bernard convection and Taylor-Couette flow. Its name comes from the field of 
superconductivity, where it models phase transitions of materials between superconducting and 
non-superconducting phases. 

In [33], Fibich and Levy showed that as €±,€2,^3 — > 0, the collapse dynamics is governed, to 
leading order, by the reduced equations (65) with 

Therefore, Continuation Results 4 and 5 hold also for the CGL continuation of the critical NLS. 

11. Continuation in the linear Schrodinger equation 

It is well known that in the linear Schrodinger equation, under the geometrical-optics approxi- 
mation, a focused input beam becomes singular at the focal point. When, however, diffraction is 
not neglected, the focused beam does not collapse to a point, but rather narrows down to a positive 
diffraction-limited width, and then spreads out with further propagation. Therefore, diffraction can 
play the role of "viscosity" in the continuation of singular geometrical-optics linear solutions. In 
what follows, we compare this continuation with those in the nonlinear case. 
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Consider the d-dimensional linear Schrodinger equation 

2iA:oV't(x, i) + AV' = 0, x = (xi, x^), (109a) 

with a focused Gaussian initial condition 

V;o(x) = e--V2g-*feorV2F^ ^ ^ 1^1^ ^-^QQ^^ 

where F > is the focal point. We can look for a solution of the form 

^(x,t) = A(x,t)e^*=o^('''*\ (110) 
where A and S are real. Substitution in equation (109) gives 

{VSf + 2St-^^^ = 0, (111a) 

and 

{A'^)^ + VS -V {A^) +A^AS = 0, (111b) 

with initial conditions 

5(x,0) = -|^, A\^,0) = e-^\ (112) 

Since ko ^ 1, we can apply the geometrical-optics approximation, and neglect the diffraction 
term A A. In this case, equation (111a) becomes 

{VSf + 2St = 0, (113) 

while equation (111b) remains unchanged. 

The solution of equations (111b, 113), subject to the initial conditions (112), is given by 

Therefore, under the geometrical-optics approximation, the solution of (109) is given by 

Cn(i'X) = ^^e"^^e''°'^vJi , 0<t<F, (114a) 

Lgo (t) 

where 

Since liuit^p P = llV'olli " '^(x)) geometrical-optics solution becomes singular at the focal 
point t = F. 

Equation (109) can also be solved exactly, (i.e., without making the geometrical-optics approx- 
imation), yielding 

i^Ut,^) = -^^e"^^e*^°^^+*^(*^^o), (115) 
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where 

L{t;h) = J- ^^ T{t;k 



1 (t tmin)'^ , t2 /j. i ^ ^ 

p — + ^Ln. ^(i;^o) = -^ 



atan ( ■ 1 + atan I — 



(116) 



and 



tmin — ^ 772 /i.2 ' -^min 



Since L{t;k()) does not shrink to zero at any t > 0, ipun exists for ah < t < oo, and in particular 
for t > F. 

It is easy to verify that the hmiting width of ipun is given by 

hm L{t;ko) = \Lgo{t)\. (117) 

ko—^oo 

In addition, since linifco^oo tmin = F and hnifco^oo Lmin = 0, then 

hm T{fM) = { \ ^ IVfl^^ (11^) 

fco->oo l~2 r < t < CO. 

Therefore, 
Lemma 14. 

f ijf° (i, x) 0<t< F, 

hm^,,„(t,x) = ^-:i/ p",^ (119) 
ko^oo ^ [Viin) (2-F-i,x)e '2'^ F<t< 00. 

Hence, hmfcg_^oo "0 coincides with tpf"^ before the focal point. Beyond the focal point, there is a 
bounded jump in the limiting phase, and the continuation is symmetric with respect to = F. 
This symmetry is to be expected, since the linear continuation is invariant under (35). 

By (118), the limiting phase is unique, both before and after the singularity. This is the opposite 
from the nonlinear case, where the limiting phase is non-unique beyond the singularity. We thus 
conclude that 

Conclusion 2. The post-collapse non-uniqueness of the phase is a nonlinear phenomena. 
11.1. Simulations 

In order to illustrate these results numerically, in Figure 23(a) we plot L{t\ ko), and observe that 
it approaches \Lgo{t)\ as fco 00, both before and after the singularity point at t = F. In Fig- 
ure 23(b) we plot r(t; ko) = avg^p{t, 0; k^) as a function of t, and observe that as ko increases, r(t; /cq) 
approaches the step function (118). 



12. Discussion 



In this study, we presented several "old" continuations and four novel continuations of NLS 
solutions. At this stage, it is not clear whether any of these continuations is the "physical" one. In 
fact, it is possible that there is no "universal continuation", i.e, that different physical setups calls 
for different continuations. 

This study suggests that all continuations share the property that the post-collapse phase be- 
comes non-unique. Indeed, this property follows from the fact that the phase of the singular 
solution becomes infinite at the singularity, and is thus independent of the specific continuation 
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2 

Figure 23: Solution of the linear Schrodinger (109) with d = 1 and F = 4. (a): Width L{t), see (116), for ko = 10 
(dashes-dots), kg = 20 (dots) and ko — 35 (dashes). Solid line is \Lgo{t)\ — |l — -Irl. (b): arg i/)(t, 0; fco) as a function 
of t for ko = 50 (dashes), fco ~ 100 (dots) and ko — 500 (dashes-dots). Solid line is (118). 



which is used. Therefore, even without knowing the "correct" continuation, we can conclude that 
interactions between post-cohapse filaments are chaotic. 

Since the Bourgain-Wang blowup solutions are unstable, they are typically classified as "non- 
generic" solutions. This study shows, however, that these solutions are "generic", in the sense that 
they arise as the sub-threshold power continuation of NLS solutions, both before and after the 
singularity. 

It is instructive to compare this study with [4]. In [4], Merle, Raphael and Szeftel showed that 
the Bourgain-Wang blowup solutions lie on the boundary of an open set of global solutions that 
scatter forward and backward in time, and also on the boundary of an open set of solutions 
that undergo a loglog blowup in finite time. This result follows immediately from the reduced 
equations (26)-(27) of the sub-threshold power continuation, see Section 4.2, since 

• The Bourgain-Wang blowup solutions correspond to /3(0) = 0. 

• For any /3(0) < 0, L(t) remains strictly positive and satisfies limt_j.-|-oo L{t) = 0. 

• For any /3(0) > 0, L{t) goes to zero at some finite Tc, at the loglog rate. 

An interesting difference between the approaches of this study and [4], is that Merle, Raphael and 
Szeftel start from the Bourgain-Wang solution at the singularity time Tc, and then find a smooth 
deformation such that the deformed solutions belong to the above two open sets on either side of 
the Bourgain-Wang solution. In this study, we start from a generic initial condition tpQ = K ■ -F(x), 
and obtain the Bourgain-Wang solution as K — > Kt^.^ 

The 'global" picture that emerges from this study is as follows. Consider a stable singular solu- 
tion of the critical NLS that undergoes a loglog collapse. Since the singular core ipji(o) approaches 
a (5-function with power Pcd it can be continued with a (^-function filament with power Per- Such a 
filament singularity can occur when the collapse-arresting mechanism leads to focusing-defocusing 
oscillations. Since this is the generic effect of conservative perturbations of the critical NLS (such 
as nonlinear saturation), see [9, Section 4.1.2], we expect that continuations that are based on 
conservative perturbations of the NLS will lead to a filament singularity. 

When the collapse-arresting mechanism is non-conservative (e.g., nonlinear damping), the so- 
lution defocuses (scatters) after its collapse has been arrested, since its power gets below Per- 
Therefore, we expect that continuations that are based on non-conservative perturbation of the 



*Of course, another difference is that, unlike this study, the results of [4] are rigorous. 
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NLS will lead to a point singularity. In addition, the same arguments as in the proof of Continua- 
tion Result 2 suggest that non- conservative continuations of solutions that undergo a loglog collapse 
have an infinite-velocity expanding core. 

In Continuation Result 2 we saw that when the continuation leads to a point singularity and is 
time-reversible, the continuation is symmetric with respect to Tc (Property 1). For this to occur, 
however, the continuation should be conservative (in order to be time reversible), yet it should 
not lead to focusing-defocusing oscillations. While this holds for the sub-threshold power and 
the shrinking-hole continuations, it is not expected to hold for conservative perturbations of the 
NLS, which generically lead to a focusing-defocusing oscillations, hence to a filament singularity. 
Therefore, we expect that continuations which are based on perturbations of the NLS equation are 
asymmetric with respect to Tc. Hence, we believe that Property 1 is non-generic. 
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